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Scalar-tensor theories of gravitation have recently regained a great interest after the discovery of 
the Chameleon mechanism and of the Galileon models. The former allows, in principle, to reconcile 
the presence of cosmological scalar fields with the constraints from experiments at the Solar System 
scale. The latter open up the possibility of building inflationary models that, among other things, do 
not need ad hoc potentials. Further generalizations have finally led to the most general tensor-scalar 
theory, recently dubbed the "Fab Four" , with only first and second order derivatives of the fields in 
the equations of motion and that self-tune to a vanishing cosmological constant. This model has a 
very rich phenomenology that needs to be explored and confronted with experimental data in order 
to constrain a very large parameter space. In this paper, we present some results regarding a subset 
of the theory named "John" , which corresponds to a non-minimal derivative coupling between the 
scalar field and the Einstein tensor in the action. We show that this coupling gives rise to an 
inflationary model with very unnatural initial conditions. Thus, we include a non-minimal, but 
non-derivative, coupling between scalar field and Ricci scalar, a term named "George" in the Fab 
Four terminology. In this way, we find a more sensible inflationary model, and, by performing a 
post-newtonian expansion of spherically symmetric solutions, we derive the set of equations that 
constrain the parameter space with data from experiments in the solar system. 
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I. INTRODUCTION 

The Galileon theory has recently emerged as an effec- 
tive theoretical realization of the Dvali-Gabadaze-Porrati 
model (DGP) 1 . The subsequent developments eventu- 
ally led to the definition of the most general second or- 
der scalar-tensor theory that includes, besides the usual 
terms of scalar-tensor or f{R) theories (as in, e.g. [31 [S]), 
also non-minimal derivative couplings to the curvature 
[TT|, 112] . This general theory, explored for the first time 
by the pioneering work of Horndeski many years ago |14| , 
provides a wide framework that virtually encompasses all 
scalar models analyzed so far in the literature. 

One interesting aspect of the non-linearities of the 
scalar sector is that they trigger the Vainshtein mech- 
anism that makes possible to build viable cosmological 
models with sufficiently small effects at local scales to 
evade Solar System constraints. Therefore, we have at 
hand an interesting alternative to chameleons [15j , in 
which the parameter space for local gravity does not 
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overlap with the one allowed by cosmic acceleration, as 
recently shown in 16J. In addition, inflationary phases 
are permitted without the introduction of ad hoc scalar 
potentials, making these models more "natural". Non- 
linearities are also responsible for new phenomena in the 
dark sector including, for instance, sub/superluminality 
and/or effective violation of the Null Energy Condition, 
thus allowing for a stable and well-defined phantom-like 
phase [ISIITS]. 

Generally speaking, the Galileon model has opened the 
way to new models for cosmology, including inflationary 
or late-accelerated ones [2TI - [23] . In fact, almost all sort of 
cosmological scenarios are possible depending on which 
Galileon model is chosen, see e.g. pS] . 

In this paper, we wish to focus on a subclass of models 
dubbed "purely kinetic gravity" [25l - f28] . or also "John" in 
the "Fab Four's" terminology [55]. This also appears as a 
special case of the very general Galileon class considered 
in M- 

We begin with a review of the equations of motion in 
a flat and empty Universe. We then complete the anal- 
ysis made in |26j by computing the number of e-folds in 
Sec. |II A| We show that a kinetically driven inflationary 
phase requires highly transplanckian values for the ini- 
tial field velocity, which rule out the model. Moreover, 
in Sec. |IIB| wc provide for a detailed analysis of the no- 
ghost and causality conditions during cosmic evolution, 
and we show that the theory becomes acausal for such 



transplanckian values unless the coupling constant is van- 
ishingly small or if the initial Hubble constant is trans- 
planckian. In passing, we prove that some claims made 
in the literature regarding the sign of the coupling con- 
stant are wrong. We argue that this model is to be dis- 
carded also because there are no reasons why the scalar 
field should be generated at all in the first place. Indeed, 
although (y9 = is only one possible solution of the equa- 
tions of motion in the vacuum, the theory is forced to 
reduce to pure GR when the same equations are solved 
inside a compact body. 

In Sec. Ill, we extend our considerations to a more gen- 
eral model, in which we include a coupling of the scalar 
field to the Ricci scalar (and thus, eventually, to matter). 
This corresponds to the "John plus George" combination 
in the Fab Four terminology. The equations of motion, 
together with no-ghosts and causality conditions, are de- 
rived and studied numerically for a cosmological back- 
ground. We find that infiation is possible provided both 
coupling constants are positive. In Sec. IIIB we derive 
the field equations in spherical and static symmetry, and 
we put some constraints on the parameter space using So- 
lar System tests. We conclude in Section IV with some 
remarks and open problems. 



II. JOHN LAGRANGIAN 

We begin our analysis by considering the action 
R 
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where R is the curvature scalar, C^^ the Einstein ten- 
sor, -0 collectively denotes the matter degrees of freedom 
coupled to the metric g^^, and k = SttGn- In this paper 
we use the mostly plus signature; 7 is a dimensionless 
parameter whereas tp has the dimension of an inverse 
length. This action is a special case of the generalized 
Galileon one presented in |53], as can one can see by set- 
ting K{X) = X, G3 = 0, G4 = 1/(2k), G5 = Kjip/2. 



A. Inflation with John 

As it was realized in [55], this model allows for a quasi 
de Sitter infiation with a graceful exit without the need 
for any specific scalar potential. Inflation is essentially 
driven by kinematics and it crucially depends on the ini- 
tial high velocity of the field, as we will shortly see. Al- 
though, in principle, the infiationary solutions begin at 
t = —00 (see [26j), we will consider the action as an ef- 
fective model only valid from few Planck times after an 
unknown transplanckian phase. Our first concern is to 
establish whether the model accommodates an inflation- 
ary phase together with reasonable assumptions for the 
initial conditions at that time. This section thus com- 
plete the analysis found in f26j by providing the number 



of e-foldings as a function of the free parameters of the 
theory. The cosmological equations in vacuum derived 
from Eq. (IT]) read 
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(2b) 
(2c) 



where a is the scale factor and a = Ina. The system can 
be partially decoupled to allow for a numerical integra- 
tion whose results are shown in Sec. II C| Isolating the 
second order derivatives, we find 
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which can be integrated numerically in a straightforward 
way. The effective equation of state (EoS) for the scalar 
field can be obtained from its stress-energy tensor or, 
more simply, by comparing our equations of motion di- 
rectly to the standard Fricdmann equations. After some 
algebra we find 



(2 + 3^27^2) (1 - ^27(^2^ 
2 + 3k27(^2 ^ 3^472(^4 ' 
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which is plotted in Fig. IT] for both positive and negative 
7. For both signs of 7, the EoS tends to —1 in the high 
energy limit (k|(/3| ^ 1), so that a large initial velocity 
for the scalar field will result in a quasi de Sitter phase. 
However, only the case of positive 7 can lead to infiation. 



Indeed, Eq. (2a) can be inverted to 3d = nip 1(2 



3k 793 ), which needs to be positive since the Hubble 
constant is a real number. Thus, 7 < implies that 
\tp\ < •^/— 2/37, which, in turn, means that we always 
have w^ > 0. Therefore, the scalar field cannot even 
start in the w^p < region if 7 < 0. More generally 
speaking, accelerated phases driven by a scalar field in 
this model require 7 > 0. 

In view of these considerations, in the following we 
shall assume that 7 is positive. Eq. (3b), t ogeth er with 



the condition that 1 — 3K7d > (see Sec. |IIB[ ), shows 
that (/3 < for an initially expanding Universe (d > 0). 
Hence, the velocity of the field decreases with time and 
uj^p is driven towards lo — +1. We characterize the end 
of inflation by the instant iendj at which uj^ = — l/Srl 



^ Throughout this paper we assume a vanishing cosmological con- 
stant. 
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Figure 1. Loglinear plot of oj^ as a function the dimensionless 
variable Kip for 7 = 1 and 7 = — 1 (dashed). 



We finally impose that inflation lasts for a number of 
e-folds N = ln(acnd/ai) greater than 60. This gives a 
relation between the initial velocity of the field (pi and 7, 
namely 



^.'>%'exp(360), 

K^7 



(12) 



which is the crucial condition for a successful (purely 
kinetic-driven) inflationary phase. We see that it involves 
a rather unusual very large pure number. In order to 
discuss naturalness, Eq. ^ is also of interest, as it fixes 
the Hubble constant at the beginning of the inflationary 
phase Hi ~ l/Zy/nrj. Therefore, the last equation might 
also be written as 
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Under the assumption that K^yip^ is initially large, one 
may derive an analytical (approximate) solution for the 
scale factor a{t) = e"'-*' and the scalar fleld at early times. 
First, we write Eq. (2a) as 3d^ = k(^^/(2 + Sn'^^ip^) so 
that 



H = a 



1 
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(5) 



where H is the Hubble constant. Integration yields the 
approximate scale factor 



a ^ tti exp 



t-U 



(6) 



Using Eq. (Is]) in Eq. (3b I, and expanding according to 
n'^ip'^ ^ 1, gives ip/ip ~ —I/^kt/ and 



ip '^ ifi exp 



t-U 

Y^K7 
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Now recall that inflation ends at tend such that 

W^(tend) = —1/3. This COlTeSpOuds to n'^Jip'^ = Cond, 

where 
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Ce„d= g (3 + V57) «1.76, 
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as can be shown by solving Eq. (B. Using Eq. (l7|, the 
condition n'^jip'^ — ^ond reduces to 



K^jipf exp 
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end; 



from which one finds 
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Replacing this in the expression ^ for the scale factor 
leads to 



'^cnd 



t,end 



(11) 



It follows that the "natural" initial conditions Hi — 0(1) 
and ifi = 0(1) in Planckian units are not allowed. On 
the contrary, a natural value for the initial expansion 
Hi — 1 (and thus 7 ~ 0.11) requires an extremely high 
transplanckian value for the initial velocity of the field 
ifii ^ 10^^ in natural units. 

It is not even possible to obtain a Planckian value for 
the initial velocity in this model, since, in any event, the 
initial Hubble constant will be greater than the one today. 
This implies that in such an inflationary scenario, ^/kt/ 
must be less than the Hubble radius today, still implying 
a very unnatural bound for the initial velocity, namely 
ipi > 10^^ in natural units. 



B. Theoretical constraints 

The characteristic feature of Galileon models is the 
derivative coupling of the scalar fleld to the metric. This 
implies a direct coupling between scalar fleld and metric 
degrees of freedom, or, in other words, that the scalar 
fleld propagation explicitly depends on the metric back- 
ground and vice- versa. Therefore, there might exist back- 
grounds for which the propagation becomes pathological 
(non hyperbolic, ie. non causal, or carrying negative en- 
ergy, ie. ghosts). In the following, we restrict ourselves 
to the cosmological background in a flat universe with 
line element 



ds^ = -dt^ 



^(*)dx2 



(14) 



and we explore the conditions for the theory to be well- 
deflned, for both scalar fleld and metric perturbations. 
We start with the scalar fleld, whose action can be writ- 
ten as 



^ip — 



a^dtd^xQ^ ip"^ 
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where 
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as it needs to be positive for the scalar field to carry 
positive energy. The field propagates at a squared speed 
given by 



1 



K7 (2a + 3d2) 
1 — Skjo^ 



>0. 



(17) 



The condition c^ > is necessary to ensure that the 
scalar field equation of motion remains hyperbolic, ul- 
timately expressing its causal behavior regardless of 
whether the scalar field perturbations are sub or super- 
luminal [171 ISOl- The two conditions above are indeed 
equivalent to the requirement that the effective metric 



n^^ 
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lif 



n^G^"" 



-ry^.'P = 0, 



(18) 
(19) 



along which the scalar field propagates, is hyperbolic with 
the same (mostly +) signature than 17^,^. 



The conditions ( 16 ) and ( 17 1 are best analyzed in terms 
of the reduced dimensionless variables x{t) = Kip and 
y{t) = ^/kq.. The first one requires 72/^ < 1/3 and it 
is automatically satisfied if 7 < 0, but also if 7 > 
provided that there is a maximal value for the Hubble 
constant H = a < 1/^/8^7. The second one reduces to 
an algebraic condition with the help of Eqs. (pi) and (Is]), 
namely 



1 - 37y' 



1 - 972/2 + 54722/-^ 



> 0^ 



(1 + 7a;2) (2 + 37a;2 
2 + 37^2 + 372a;4 



>0. 



(20) 

In summary, Eqs. (16 17) require 7 > and \y\ < 1/y/S^ 

< 



or 7 < and the two possibilities |a;| > l/\/— 7 or |a;| 
•\/— 2/37. These conditions are therefore less restrictive 
than the one implied by the Friedmann equation (see 
previous section). 

The curvature background implies non-standard prop- 
agation for the scalar degree of freedom. In a quite simi- 
lar fashion, the scalar field background modifies the stan- 
dard spectrum of metric perturbations. Similar condi- 
tions for the avoidance of ghosts and euclidean metrics 
also exist. These have been derived in full generality in a 
very wide class of Galileon models in [3l] , whose conven- 
tions we follow. These conditions, namely Eqs. (23), (25- 
27) of [3T] reduce to rather simple algebraic constraints 
in our case, after the necessary manipulation using the 
equations of motion ^ and (|3| 



Qt > 4^ 1 + ^ > 0, 

2 

I - 2 - 



(21a) 
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the scalar part of the metric perturbations, namely 
Qs>0-^ TT-r^ — ^ > 0, 



c| > 0^ 



2 -I- 37x2 

12 + 36jx^ + 19-f'^x'^ - U'j^x^ - 37^x8 
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This whole set of equations is difficult to reduce alge- 
braically because of the last one. However, one might 
easily plot the six functions of x defined above, and one 
typically finds that both positive and negative values for 
7 are allowed on a given range |x| < ^^, where typic_ally 

^-y behaves as O 



(l/\/|7l)>seee.g. Figs. 2 and Fig. 0be- 



low. Hence, large (transplanckian) valuesior |x| are only 
allowed for small I7I ^ 1. This means that the space for 
possible velocities of the field x = Kip needs to be typ- 
ically subplanckian, unless 7 is vanishingly small. This 
will be linked to the results found earlier, where trans- 
planckian initial velocity were required for a successful 
inflation, leading to negative squared speeds c| and c^ 
in that epoch. This is shown in Fig. [5] 




Figure 2. Plot of the six conditions Qi and d derived above, as 
a function x. Here 7=1. Allowed values for the field velocity 
is typically |a;| < 0(7"^"), as made clear by drawing similar 
plots while varying 7. 




where Qt and ct are defined as their scalar counterpart 
Q^p and c2 above, but stand for the tensor perturbations 
of the metric field. Similar conditions need to hold for 



Figure 3. Plot of the six conditions Qi and Ci derived above, 
as a function x. Here 7 = — 1. Allowed values for the field 
velocity is typically \x\ < 0(171"^/^). 



In passing, we note that the claim made in the hter- 
ature (see e.g. [24l [281132]) according to which only the 
subclass 7 < is a ghost-free theory is wrong (at least in 
the background considered here). Notice that the scalar 
field is well-defined although being a phantom in certain 
regime (in the case 7 < 0), a situation reminiscent of the 
one discussed in [f 8'. However, as shown previously, the 
Friedmann equation actually prevents the scalar field to 
enter this regime. 



C. Numerical results 

In this section, we quickly show the cosmological behav- 
ior in the John model, for both positive and negative 7. 
As discussed before, the negative 7 case leads only to 
a decelerating Universe: the phantom regime is not an 
acceptable initial condition (as it entails an imaginary 
Hubble constant), and neither can be reached. Only pos- 
itive 7 leads to acceleration, and to an inflationary phase 
in the early Universe, a drawback being the presence of 
non-causal behavior for the scalar and tensor perturba- 
tions of the metric. These plots (Figs. |4]-[7| have been 
obtained by numerical integration for an initial condi- 
tion of Lpi — 10 in natural units in the case 7 = 1, and 
ifi = 0.1 for 7 = —1. 
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Figure 5. Evolution of the parameters for metric perturba- 
tions (scalar and tensorial parts), respectively Qs, c|, Qt, 
Cy, as defined in the text vs the scale factor, for 7=1. As 
found theoretically, the initial high velocity of the field drives 
both the speed of (scalar and tensor) metric perturbations to 
imaginary values, thus signahng a breakdown of hyperbolicity 
for metric perturbations. (Here it happens when the corre- 
sponding curves terminate, since the y axis is in logarithmic 
scale). 
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Figure 4. Evolution of the acceleration parameter q = -\-aa/a?' 
vs. the scale factor, for 7=1. Also shown in the graph: the 
evolution of the scalar field EoS m.^ and Q^ and cj. The 
scalar field becomes superluminal during the transition from 
a de Sitter Universe to a to stiff matter-dominated one. How- 
ever, as discussed in the text, it remains hyperbolic and, thus, 
causal. Moreover, it carries positive energy. 



D. Discussion 

We have established that kinetically driven inflation 
in the Galileon theory involving the simplest coupling 
to Einstein tensor is not viable. It requires unnatural 



Figure 6. Evolution of the acceleration parameter q vs. the 
scale factor, for 7 = —1 with initial condition satisfying \xi\ = 
The field starts with an EoS w^ 



0.1 < \/-2/37. 

the Universe only decelerates. 



0.5 and 



transplanckian values for the initial velocity of the field, 
which, in turn, implies various instabilities, as discussed 
above. 

This model has anyway another serious drawback. In 
the absence of any direct coupling to the Ricci scalar, 
there is no reason why the scalar field should be generated 
at all (even in presence of cosmological matter fluid) . In 
other words, (/3 = is always a solution in this class of 
models, whatever the matter content is. At local scales. 
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Figure 7. Evolution of the parameters for metric perturba- 
tions (scalar and tensorial parts), respectively Qs, c|, Qt, Ct, 
as defined in the text vs the scale factor, for 7 = — 1. This 
model with 7 < is well-behaved but do not accommodate 
inflation. 



this problem appears in the following way. We checked 
numerically that a relativistic star (with flat asymptotic 
conditions) must be described by pure GR, i.e. ip = 
for all r is the only solution that is regular at the center 
of the star. To conclude, the model considered so far is 
trivial in the sense that it cannot be different than GR, 
except if one imposes non- vanishing initial conditions for 
the scalar field at early times. 

This acts as a leitmotiv for a more realistic model, 
studied in the next section, where we add to the La- 
grangian a direct coupling to the Ricci scalar: R{1 + £1^9), 
namely the " George" term in the Fab Four terminology. 



III. GEORGE AND JOHN 



We now consider the extended model given by 
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where e is a dimensionless, free parameter. Of course this 
is not the most general coupling one might consider but 
it is anyway reminiscent of Brans-Dicke coupling. Notice 
that one might worry that the effective gravitational con- 
stant Geff = G/(l -I- €^/k,lp) might easily become negative 
in this model, meaning that the action chosen here shall 
trivially lead to dynamical pathologies for cLp sufficiently 
large and negativq^ 



Such an argument would call in favor of defining a 
better coupling function F in the Georges term F{(f)R. 
However, this would be a misleading conclusion here, 
since the John term couples the derivatives of the metric 
and of the scalar field, thus impacting their propagation. 
Therefore, only the entire set of no-ghost conditions to- 
gether with causal propagation conditions (positivity of 
the squared velocities) for both the scalar and the metric 
perturbations can decide which regions of the configura- 
tion space are well-behaved. This is done in the following 
sections (on a cosmological background), based on the 
conditions derived in AppendixlB] In this light, the func- 
tion F chosen above is just the simplest one could have 
chosen, and might furthermore be understood as retain- 
ing only the first term in a series expansion of a more 
general function F. 

The cosmological evolution in this theory is typically a 
function of four parameters, the initial value of the field, 
its velocity, and of the two dimensionless parameters 7 
and e. It goes beyond the scope of the present paper 
to provide with a comprehensive study of this parameter 
space. However we highlight some essential features of 
the model thanks to numerical results displayed in the 
next section. We provide both the cosmological evolu- 
tion and the analysis of causality and positivity of energy 
within subclasses of the model, depending on the signs 
of 7 and e. 



A. Cosmological behavior 

The equations of motion in a flat, empty Universe, 
derived from Eq. ( 23 ) are given in Appendix Kj see 
Eqs. (A1||A3). We extended the analysis of the no-ghost 



and causality conditions to this more general framework, 
and we also provide the scalar field EoS, see Appendix [A| 
andE 

The numerical results are the following. The case e — 1 
and 7 = 1 is pretty similar to the case John alone, see 
Figs. (Is]) and ^ below. Inflation thus occurs in the case 
e > and 7 > 0, but the acausal behavior still shows up 
in the very early Universe. The number of e- folds is a 
function of the two initial conditions for the field and its 
velocity, and the dimensionless parameters e and 7. A 
further analysis that goes beyond te scope of this paper 
would determine whether the addition of the George term 
help in solving the naturalness problem encountered with 
John alone in Sec. HIl 

The case e = — l,7 = lis clearly pathological for the 
various no-ghosts and no-acausal conditions, as seen on 



^ In fact, what matters in the case 7 = 0, is that the scalar field 
propagates positive energy in the Einstein frame. Performing a 



conformal transformation, this is equivalent to the usual Brans- 
Dicke condition 2ii) -|- 3 > 0, where w = e^(l + e^JHip) here. 
Then, our model with 7 = would indeed be pathological if 
if> < — (3/(2e^) -I- 1)/{€^/k) . However the 7 term introduces 
new terms in the equation of motion for the scalar field which 
invalidate such a conclusion in the general case 7 5^ 0. 
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Figure 8. Evolution of g, vj^, Q^ and cj with the scale factor. 
We observe the same transition from inflation to stiff matter 
for the scalar field. Initial conditions are ipi — fOO and ipi = 1. 
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Figure 9. Evolution of Qt, Qv and c^jC^ with the scale 
factor. These two speeds are negative in the early universe. 



Fig. (10). Actually this theory leads to a double infla- 



tion scenario (see the acceleration parameter): the Uni- 
verse transits from one de Sitter phase to another one, 
and experiences in between a super-acceleration phase. 
Finally, the case with negative 7 is similar to what we 
found for John alone: the theory is well-defined, ghost 
free and causal, but fail to exhibit any acceleration at 
all, see Figs. 



(10) 



George" . 

In a first step we solve asymptotically the field equa- 
tions in spherical symmetric and isotropic coordinates. 
The line element reads 

ds^ = -~A{r)^dt^ + B{rY (dr^ + r^d^D) , (24) 

and the equations of motion for the components and for 
If are given in Appendix [Cj The post-Newtonian analysis 
begins with the expansion of the metric components as 
A^ = 1 + E^a^AN 52 = 1 + Y.^b^/r' and if = Po + 
'YliiPil'^^ and of the field equations in powers of 1/r. By 
equating the coefficients of equal powers of r we find 
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where r, = 2'^.*^ is the Schwarzschild radius of the cen- 
tral body, z — 1 + ^, and 7 = 7^2. In the expan- 
sion above, we neglected higher order terms in rg/r^ in e, 
in pi/{rMp) and in ^/{rMp) (which means we suppose 
these terms to be smaller than 1). We recall that, in our 
conventions, 7 and e are dimensionless parameters. The 
asymptotic scalar field value po = fi''' -^ 00) (in GeV) is 
a free parameter that can eventually be connected to the 
cosmological evolution of the scalar field. 

The next step then amounts to relate the dimensionless 
parameter pi to the scalar charge of the central body, and 
thus eventually to the couplings e and 7. This could be 
determined numerically by solving the equations in the 
interior of the body [33] with suitable asymptotic and 
regularity conditions, together with, e. g. a polytropic 
equation of state for the star's interior. Such a discussion 
about possible spontaneous scalarization goes however 
beyond the scope of the present paper and is left for 
future studies. In this section we only sketch how solar 
systeme constraints can be derived. 

The final step consists then in computing observables 
effects from the metric Eqs. ( 24p5 l, which show a de- 
viation from standard General Relativity in the Solar 
System. A first constraint comes from the anomalous 
perihelion shift. First of all, we determine the geodesic 
equation (to first order in the metric deviation) to find 
the planetary equations of motiorrl The anomalous term 
is treated perturbatively to find the perihelion shift with 
the Gauss equations, which determines the evolution 
of the different orbital parameters due to perturbative 
forces. In particular, we derive the rate of change of the 
argument of the perihelion uj. Finally, the secular term of 



B. Solar system 



In this section we show how to derive solar system 
constraints on the free parameters of the model "John -|- 



^ When computing the geodesies, we suppose that matter is min- 
imally coupled to the metric. 
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Figure 10. Cosmological evolution of q,Q^,c^,Qs,(?s,QT,CT with the scale factor. (Top); the case e — —1 and 7 = 1 
corresponds to a "double inflation". The scenario is pathological in many respects: singularity in the scalar field velocity (see 
top left) while the scalar field EoS is imaginary. Also (see top right, the y axis is logarithmic): c^ < and there are periods 
for which c^ < and Qt < 0. (Bottom): the case e — 1 and 7 = —1. The model is well behaved but does not accelerate the 
expansion, ie. as in the John alone model. The universe is actually in a super-stiff regime, and hence in a highly decelerating 
phase. 



the change of the argument of the perihelion is exhibited 
by a time average (over an orbital period T) 
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The results of this whole procedure is given by 
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4ars(l - e2) 
where a is the semi-major axis of the orbit, e its eccen- 

T 



tricity, and n its mean motion n = ^. The advance 



of perihelion of Solar System planets is very tightly con- 
strained by planetary ephemerides. In particular, IN- 
POPlOa gives constraints on supplementary advances of 
perihelia (see Table 5 of [311 )■ Table IT] gives the value of 
the tti coefficients appearing in the expression of the ad- 
vance of perihelion (27 1 and the constraints coming from 
INPOPlOa. 

The most stringent constraints in the case considered 
here are obtained by data from Mercury and read 



3.12 X 10^^ m* < Jj- < 6.25 x 10^° m^, (29a) 






-2.06 X 10^'^ m* < ^i^ < 4.12 x 10^"^ m*, (29b) 
-1.13 X 10^° m^ < j^ < 2.26 X 10^ m^, (29c) 
-5.16 X 10^2 ^ < 7^- < 1-03 X 10^2 m. (29d) 
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Figure 11. Incompressible anomalous signal due to the dif- 
ferent terms in the metric ( 24]|25 I on Cassini Doppler signal 
with indicated values. The produced residuals are larger than 
Cassini accuracy (which is 10~^*), therefore these values for 
the parameters are ruled out. 



accuracy (10 ^'^) gives the boundary values 



These constraints are obtained by considering only de- 
viations on the dynamics (i.e. on the equations of mo- 
tion). Other constraints can be derived using propaga- 
tion of light rays in the solar system. For example, radio- 
science experiments include light propagation through a 
Shapiro-like term that can be derived from the expres- 
sion of the metric ( 24p5 ). One way to obtain such con- 
straint is to follow the strategy and to use the software 
of [351 - 137] . The main idea presented in these papers is to 
simulate radioscience observables directly from the space- 
time metric so that the software includes both deviations 
on the dynamics and deviations on light propagation. In 
particular, we use this software to simulate a two-way 
Doppler link between Earth and Cassini spacecraft from 
May 2001 on the metric ( 24p5 l and to analyze them 
in GR by fitting the initial conditions of the spacecraft. 
The residuals that emerge are the incompressible signa- 
ture produced by the alternative theory considered on 
Doppler signal for Cassini. Comparing this signal to the 
Doppler accuracy of the mission (10~^^) allows us to give 
order of magnitude of constraints on the theory. Fig. [TT] 
represents the incompressible signatures produced by pa- 
rameters entering the metric ( 24 



25 1 on Cassini Doppler. 
The three sharp peaks occurred at solar conjunctions. 
The order of magnitude of the residuals observed in this 
figure is larger than Cassini accuracy, which means that 
the values of the parameters should be smaller than the 
indicated values. 

We run a set of simulations with different values for the 
parameters appearing in the metric ( 24p5 |. Fig. 12 rep- 
resents the evolution of the maximal Doppler residuals 
obtained in Cassini signal as function of metric param- 
eters. Requiring the residuals to be lower than Cassini 
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< 5.56 X 10" 



(30d) 



It should be noted that radioscience constraints are sig- 

2 

nificantly better for j^j-^ and -^^^ as compared to peri- 
helia advances. On the other hand, the constraint from 
the ephemerides on -j^r^ is significantly better while the 
constraint on j^ is of the same order of magnitude. 

We thus see that solar systems tests of GR will put 
severe constraints on the free parameters of the model. 
As already stressed however, the full numerical study of 
Solar's interior is required, while a thorough analysis of 
possible spontaneous scalarization (depending also of the 
Georges coupling) is necessary. 



IV. CONCLUSIONS 

In this paper, we have explored some phenomenology as- 
sociated to a subset of the "Fab Four" scalar-tensor the- 
ory. The philosophy behind this preliminary work is that 
we cannot forget about solar system constraints on the 
parameter space, even when we deal with inflationary so- 
lutions. Traditional models of inflation rely upon the fact 



10 





ai (mas/cy/m ) 


a2 {mas/cy/m ) 


0:3 {mas/cy/m ) 


04 {mas/cy/m.) 


as {mas/cy/m^) 


dw/dt {mas/cy) INPOPlOa 


Mercury 


3.2 X 10"^^ 


4.85 X 10"^'' 


8.83 X 10"" 


19.4 


8.2 X 10"'' 


0.4 ±0.6 


Venus 


7.69 X 10"=^'' 


2.53 X 10'^^ 


9x 10-^^ 


3.89 


1.65 X 10-* 


0.2 ±1.5 


Earth 


1.3 X 10"^* 


5.9 X 10"^'' 


2.9 X 10"'^ 


1.73 


7.33 X 10-^ 


-0.2 ±0.9 


Mars 


1.36 X 10"^^ 


9.13 X 10"^^ 


6.77 X 10"" 


6.1 X 10"' 


2.58 X 10"^ 


-0.04 ±0.15 


Jupiter 


1.51 X 10"=^* 


3.56 X 10"^'^ 


9.06 X 10-^'' 


2.81 X 10-^ 


1.19 X 10-*^ 


-41 ± 42 


Saturn 


5.33 X 10-''° 


2.3 X 10-^° 


1.08 X IQ-'^ 


6.16 X 10"^ 


2.61 X 10"' 


0.15 ±0.65 



Table I. Values of the diflferent ai coefficients (28 1 for planets and constraints given by INPOPlOa for supplementary advances 
of perihelia (values given in |34] ) . 
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Figure 12. Representation of the maximal Doppler resid- 
uals signal due to John±George modification of gravity 
parametrized by the four parameters in the metric ( 24][25 l 
for the Cassini mission between Jupiter and Saturn. The red 
(dashed) lines represent the assumed Cassini accuracy. 



that the inflaton field decays, at some stage, into ordi- 
nary matter through some reheating mechanism. There- 
fore, the scalar-tensor natm'e of inflationary gravity is lost 
very soon in the evolution of the Universe. On the oppo- 
site, in the models studied in this paper the scalar field 
should live and shovsr its effects until nowadays. There- 
fore, the parameter space determined by constraints from 
cosmological observations must overlap with the one de- 
termined by solar system tests. The "Fab Four" theory 
has many parameters with a very rich phenomenology, 
and with this paper we begin an ambitious plan for its 
systematic study. 

In this work, we contained ourselves to the cases John 
and John plus George. The John case represents a theory 
with a non-minimal derivative coupling between scalar 
field and Einstein tensor. It was already known that this 
model admits inflationary solutions with a graceful exit. 
Here, we constrained the coupling constant 7 by showing 
that it must be positive in order to have successful infla- 



tion. If this is the case, however, very unnatural initial 
conditions are required. In particular the field velocity, 
which is related to the energy density, must be huge com- 
pared to the Planck scale, rendering the theory no longer 
trustworthy. Negative values for 7 are permitted but do 
not allow for inflation while ghost states might appear. 
The most serious problem however comes from the fact 
that the model turns out to be trivial when one tries 
to solve the equations of motion inside a compact object. 
Indeed, we found that the only solution with a finite field 
at the centre is (/3 = everywhere. 

These facts have convinced us to extend the theory 
to include the term named "George", which is nothing 
but a coupling between the scalar field and the Ricci 
scalar. The parameter space is now two-dimensional and 
we solved numerically the equations of motion for a cos- 
mological background. The main result is that the sign 
of the two coupling constants must be positive in order to 
have an infiationary phase with graceful exit. We reserve 
the analysis of the naturalness problem, i.e. the need 
for extreme initial conditions, for future work. We also 
performed a post-Newtonian analysis of the theory by 
solving the equations of motion by imposing static and 
spherical symmetry and expanding the fields for large ra- 
dial distance. The aim was to put some constraint on the 
free parameters, which now include also the asymptotic 
value of the scalar field and the scalar charge, which now 
makes sense as there are non-trivial solution also for the 
interior of a compact object. The results still allow for a 
large parameter space, therefore future work is necessary 
in order to improve the constraints. 

As mentioned above, this paper is the first of a se- 
ries that aim at a systematic study of the "Fab Four" 
phenomenology, in order to isolated the experimentally 
viable sectors of the theory. Besides these aspects, there 
are several issues that deserve further investigations. For 
example, if these models truly lead to infiationary cos- 
mology and graceful exit, we will need an alternative re- 
heating mechanism. Maybe, the complexity of the theory 
acts as an effective potential for the scalar field, which 
resembles the usual power-law forms of usual inflation. 
But if this is not the case, then we need to find alter- 
native explanations. Another aspect that needs to be 
studied is the relevance of these modifications of gravity 
in terms of late-time cosmology, as we expect modifi- 
cations driven by the scalar field. This would include. 
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at the background level, the study of tracking solutions 
and of the convergence mechanism towards GR, if any. 
This has not yet been adressed in the litterature for the 
John Lagrangian (but see [35] )• The study of cosmolog- 
ical perturbations, in particular CMB spectra and large 
scale structures, might then further reduce the parameter 
space. Finally, gravitational effects might be relevant at 
galactic scales and give rise to alternative explanations 
to the anomalous galactic rotation curves. 

In conclusion, we believe that the recent developments 
in scalar-tensor theories of gravity have opened the door 
to new and intriguing research directions, and we are 



confident that many interesting results will be obtained 
in the near future. 
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Appendix A: Cosmological equations 

In terms of the reduced variables x{t) = Kip, y[t) = ^JKa and z{t) = 1 + e^/Kip{t), the equations of motion for a flat 
and empty universe derived from action Eq. ( 23 1 are 



6exy + x^ (-1 + Q^y^) + &y^z = 

Ax (e + 'fy/Kx) y + x'^{l + S'fy^ + 2-fy/Ky) + %y^z + 2^/ll{fix + 2yz) = 

3y {x — 2ey — S'yxy ) + y/n (i — 3jxy — 3(e -I- 2"fxy)y) = 

which can be decoupled in the following way: 

—3a; [ex + 4 (e^ -t- ^x"^) y + 7e^xy'^\ -+ 6y{—2x + ey)z 
^~ 2k1/2 (3g2 _^ i2e-ixy + 7^2 (1 + 9-fy^) + {2- 6jy^) z) 

2exy (1 - 1572/2) + x^ [~1 + 3jy^ (4 - 9jy^)] - Gy^ {2e^ + z- 3jy^z) 



y = 



The scalar field EoS is given by: 



2k1/2 [3e2 ^ Uejxy + 7^2 (1 + 9'yy'^) + (2- 672/2) z] 



= -X (37X^ + 2z) 



N 
D 



with 



N = 6x {jx^ - z) (37x2 + 2z) + eVSe^y^x^ (3e2 + 3jx^ + 2z) {l^ix"^ -f 2z) - 18e^ (77x3 -|- 2xz) 
-2\/3eV'a;2 (3e2 + 37x2 + 2z) {p.'fx^ -)- 167x2z - 4^2) -|- 9e2 (l572x^ -I- 47x'''z - 4x^2) 



D = 



-3ex + v/3 v/a;2 (3e2 -^ 37x2 + 2z) 

V&r'x^ -F ZQ-i^x^z -h 247x^^2 -h 8z3 -K 6V3e7X {jx^ + 2z) y/x^ (3e2 + 37x2 + 2z) + 3e2 {i-f'^x^ + 4^2) 

Appendix B: Ghosts conditions 



(Al) 
(A2) 
(A3) 



(A4) 
(A5) 

(A6) 



(A7) 



(A8) 



The coupling to the Ricci scalar, "George", does not change the analysis made for the scalar field sector of the theory. 
Thus the two following conditions still hold 



?^ = o (1 - 372;') > 



_ 1-7(32/^ + 20^2/) > „ 



1 - 372/2 



(Bl) 

(B2) 



For the metric perturbations, we derive, based on Eqs. (23), (25), (26), and (27) of [5T] 
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Qt>0^z+^>0 
c|>0^z-^>0 



(B3) 

(B4) 

(B5) 



for the tensorial part, and also 



„2 «„,2. 



Qs > ^ 3e^ + 12e7xy + Q^x V + 2z + j Ix^ - 6y^z] > 



(B6) 



for the scalar part of the metric perturbations, while their squared speed cj. > leads to 



2y I — h z j (ex + ijx'^y + 2yz) + 2a; (e + "f^/nx) (72;^ + 2z) (ex + 3'yx'^y + 22/z) 



1 2 

+ - (yx"^ — 2z) (ex + Sjx^y + 2yz) 



2 \ 2 r 
-2y^( ^+2 



, . 2xy . 
e X H ^ + 3"fx[2xy + xy) + 2yz 



>0 



(B7) 



Appendix C: Spherically symmetric equations of motion 



We derive the equations of motion for the action ( 23 ) with a spherically symmetric and static field configuration. We 



consider the metric (24), and we replace its components in the Lagrangian. With the Noether theorem, we find the 



equations of motion for the fields A, _B, i?, and ip. Finally, we impose the gauge R = r, and we find three equations 
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plus a Hamiltonian constraint that read 

= (2ip''^Br^B" + iip'B^rtp" + 4:ip'Br'^ip"B' - bip'^B'^r'^ + 2ip'^ B'^) -/ k^ 



- 4BV-B"r^ - 2B^B'r^Lp' - 4B^ip'r - SB-^ipB'r - 2B'^r\" 



2B 



VS'V^^e^i 



/2 



+ 



8B-^rB' - AB'^r^B" + 2B^r^B'^ - (p'^K^BV 



(CI) 



= Uip'^B'ArB + 2lp'^B^A - Sip'^B'^Ar"^ + 4(p' Br^ ip" AB' + 8ip''^B'^rA' + 2(p'^B'A'Br^ + 

+ 2ip'^A"r^B^ + 4ip'B\^"A + 4p'A'r^ip"B^ + 2ip'^ABr^B"^-fK^ + 

AB^AB'r'^ip' + 2B^AipB'^r^ - SB'^Aip'r - &iprA'B-^ - 8B^AtpB'r - 4:B^AipB"r'^ 

4ipB'A'r^B^ - GA'r^^'B'' ~ iifA"r^B^ - 4B^Ar^^"yK^/^ + 

AB^Ar'^B" - Ar^A"B^ - SrA'B^ - SB^ArB' - AA'r^B'B^ + 2B^Ar^B'^ - ^'^n^B^Ar^, 



(C2) 



== Uip'A'B^ + Aip'AB'B^ + Aip' B' Ar^ B" B + Aip' B" A' r^ B"^ + Aip"rA'B^ + Aip'rA"B^ + 
+ V'B'A'r^B^ + Stp'B'A'rB^ + 4:ip'B"ArB^ + Aif"B'ArB'^ 
- Aip'B'^ArB - Gip'B'^A'r^B - Gtp'B'^Ar^ + Aip' B' A" r'^ B^ + 2ip" B'^ Ar^ b) -/ k + 

-1/2 



- SB'ArB^ - AB^rA' - Ar^AB"B^ - 2B'A'r^B-^ - 2B''r^A" + 2B'^Ar^B^ 
2AVVB'^ + 2B^r^ip"A + Atp'B^^Ar + 2B'r^ip'AB*, 



(C3) 



= hifi'^A'B^ + 2r^'^A"B^ - 6(p'^B'^Ar + Aip' Aip" B"^ - 2^'^B'AB - Aip''^A'B'rB 

+ Aip'rB'Aip"B + Arip'A'ip"B'^ + 2ip''^rB" AB\-iK^ + (- AipB'AB^ + AripB'^AB'^ 

- AB'^rA'ifi' - AB-^ipA"r - AB^ip"Ar - AB'^ipA' ~ AB'^ip'A - A^pB" ArB^^eK^/"^ 
+ ArB''^AB^ - ArAB"B^ - AB'AB^ - AB*A' - AB^A"r - 2B^ ip'"^ Ar k^ . 
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